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ABSTRACT 

We present a comprehensive analysis of the X-ray hght curves of SN 1993 J in a 
nearby galaxy M81. This is the only supernova other than SN 1987A, which is so 
extensively followed in the X-ray bands. Here we report on SN 1993 J observations 
with the Chandra in the year 2005 and 2008, and Swift observations in 2005, 2006 
and 2008. We combined these observations with all available archival data of SN 
1993J, which includes ROSAT, ASCA, Chandra, and XMM-A/"ewton observations 
from 1993 April to 2006 August. In this paper we report the X-ray light curves 
of SN 1993 J, extending up to fifteen years, in the soft (0.3-2.4 keV), hard (2- 
8 keV) and combined (0.3-8 keV) bands. The hard and soft-band fluxes decline 
at different rates initially, but after about 5 years they both undergo a decline. 
The soft X-rays, which are initially low, start dominating after a few hundred 
days. We interpret that most of the emission below 8 keV is coming from the 
reverse shock which is radiative initially for around first 1000-2000 days and 
then turn into adiabatic shock. Our hydrodynamic simulation also confirms the 
reverse shock origin of the observed light curves. We also compare the Ha line 
luminosity of SN 1993J with its X-ray light curve and note that the Ha line 
luminosity has a fairly high fraction of the X-ray emission, indicating presence 
of clumps in the emitting plasma. 
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1. Introduction 



The explosion of a massive star as a supernova (SN) can drive powerful shocks into the 
circumstellar medium (CSM) of the progenitor. The CSM is established by the mass lost 
from the progenitor prior to explosion. Collision of the ejected material with the CSM leads 
to a blast wave shock with velocities 10, 000 - 20, 000 km/s and a hot shell of T ~ lO^K. 
The interface of the expanding ejecta where it first meets the circumstellar gas is the contact 
discontinuity which itself propagates outwards with time. As the ejecta expands, the interface 
region between the blast wave shock and the contact discontinuity sweeps up matter from 
the circumstellar gas and decelerates due to the accumulated mass. This in turn gives rise 
to a shocked shell that propagates back into the cold expanding ejecta This 'reverse' shock 
propagates inwards at ~ 10^ km/s, significantly slower speed than the fastest expanding 
stellar ejecta and has a temperature of T ~ 10^ K. The reverse shock is believed to be the 
site of most of the observable X-ray emission at late times. 

In this paper, we describe the long term light curves of the SN 1993J in various X-ray 
bands. We briefly describe in §2] the current understanding about the unusual nature of SN 
1993J and its progenitor star as obtained from multi-waveband observations and theoretical 
analysis. We discuss the X-ray emission mechanisms and the previous X-ray studies of SN 
1993J in IJHI In we discuss the X-ray observations and analysis of SN 1993J over the years 
made with several X-ray telescopes. Our detailed results and interpretation are mentioned in 
§3 Hydrodynamical simulations are discussed in §H1 and X-ray versus optical Ha correlation 
in ^ We report our main conclusions in ^ 



2. A supernova that underwent a metamorphosis of spectral types 



SN 1993J is one of the best studied supernovae (SNe) in all wavelengt h regimes, second 
only to SN 1987A. It was visually discovered on 1993 March 28. 906 UT (IRipero fc Garcia 
1993h . in the nearby galaxy M81 (aka NGC 3031 at d = 3.63 Mpc. IPreedman et al.lll994h . It 



was the optically brightest SN in the northern hemisphere since SN 1954A, having reached 
a secondary maximum brightness of V = 10.8 mag on day 21.1. 



SN 1993J is a SN that has undergone an "identity crisis". It's spectrum underwent a 
transition from a type II spectrum (characterized by strong hydrogen Balmer lines) at early 
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epochs, to a type Ib-like spectrum at ~ 300 days (with its nebular spectra having weak 
hydrogen but strong He I hues). This SN was classified as Type lib S N and provided for the 



first time a link be tween type II supernovae (SNe) and type lb SNe (IFilippenko et al.lll993 



Swartz et al.lll993l ). Models of SN 1993 J based on the early light curve indicated that the 
progenitor star lost all but a small amount of its hydrogen layers due to mass transfer in a 



binary system (INomoto et al.lll993l : iRay et allll993l : iPodsiadlowski et allll993l : IWoosley et al. 



19941 ) ■ Thus the shock heated photosphere could quickly recede through the small H layer 
into the deeper He layers during the initial expansion and cooling phase itself. Among core 
collapse SNe, it is thus possible to have a continuum of hydrogen envelope masses remaining 
on the progenitor star when it explodes. 

At early epochs SN 1993J showed th e typic al signatures of circumstellar (CS) inter- 
action in the radio (IVan Dvk et al.l Il994l ). UV (IFransson &: SonnebornI Il994l ) and X-ray 
( IZimmermann et allll994j ) wavelengths. The UV and optical spectra taken with the Hubble 
Space Telescope (HST) and Keck telescopes have revealed the signature of a massive star, 
along with the fading superno va 10 years after ex plosion, which is the binary companion of 
the progenitor that exploded ( Maund et all 2004 ). The circumstellar medium has therefore 
been modified by the mass loss in a binary system in SN 1993J. 



3. X-ray emission and early X-ray studies of SN 1993J 



The X-ray luminosity and its time variation depends, among other things, on the density 
structure of the stellar ejecta. In core collapse SNe, the explosion dynamics quickly leads to 
an ejecta outer density profile with a steep power-law ( pej r"" where n is a constant) in the 



radia l coordinate, but with a r elatively fiat inner core density profile (j Chevalier and Soker 



19891 : iMatzner fc McKedll999l ). The ejecta density profile depends on the initial structure 
of the sta r and, in particular, is aff ected by if the progenitor had a radiative or a convective 
envelope (IMatzner &: McKeelll999l ). The shock propagation through the outer profile does 
not depend upon the behaviour of the inner layers and a limiting structure is described by 
a self-similar solution as the shock front accelerates while propagating through the outer 
layers with rapidly decreasing density. The self similar nature of the evolution of the ejecta 
dominated SNe (or supernova remnants (SNRs)) applies only for ejecta with steep enve- 
l ope index n > 5, since f or such steep n, the mass and energy of the ejecta remain finite 
(iTruelove k McKeel ^9991 ) . 



The shocked shells (forward or reverse) generated due to ejecta-wind interaction have 
very high temperatures (f|T]) and can emit X-rays. The forward shock is adiabatic almost all 
the time and its luminosity follows the time dependence of ^forward cx: (I Chevalier &: Fransson 
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20031 ). However, depending upon the ejecta density profile and the mass loss rate of the SN 
progenitor star, the reverse shock can be either adiabatic or radiative, or can be radiative 
at early phase and then make a transition to adiabatic phase. The e xpression for free-free 
emissi on fr om adiabatic reverse s hock and its luminosity is derived by I Chevalier fc Fransson 
koO^ and lFransson et all Jl996h : 

(n-3)(n-4)2 



rad ^ 
rev ~ 



3 X lO'^gs 



A(n - 2) 



M_5 



t 

1 day 



ergs 



-1 



(1) 



Here M_5 is progenitor mass loss rate in units of lO~^M0yr~^, u^i is the wind velocity in 
units of 10 km s~^, and gs is the free-free Gaunt factor. This shows that free-free emission 
from adiabatic reverse shock follows time evolution of t^^ for constant mass-loss rate and wind 
velocity. H owever, if line emission doni inates, the time dependence of X-ray line luminosity 
is ii-2-2.2m (Ichevalier fc Franssonlll994j ). where m is the expansion parameter in i? oc t"*. 



Fransson et al.l (119961 ) have discussed that in the initial phase, reverse shock with high 



density gradient is likely to be radiative resulting in formation of a cooled shell between 
the reverse shock and the forward shock. When the ele ctron temperature is T ), < 2 x 10^ 
K, the line emission dominates the total X-ray emission. I Chevalier fc Fransson! (120031 ) have 
discussed the importance of line emission and its effect on the cooling rate of the gas behind 
the reverse shock. In case of line emission, a thermal instability develops and the gas cools 
up to about 10^ K where the temperature is stabilized only by photoelectric heating from 
the shock balancing the cooling. For an ejecta velocity scale of Vej, and a reverse shock 
moving through an ejecta density gradient p cx r^" with solar composition one obtains a 
cooling time scale fr om tcao] = 3kTp/nA, where A is the cooling funct ion. The cooling time 



can be expanded as (IFransson et al.lll996l : IChevalier fc Fransson! |2003| ): 

-1 
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Since large exponents of the velocity scale and density gradient index are involved, it is clear 
that the cooling time is sensitive to them, and also the mass loss rate of the pre-explosion 
progenitor star. The most important effect of cooling gas between the reverse shock and 
the observer is that the cool gas absorbs most of the emission from the reverse shock, and 
in spite of the higher intrinsic luminosity of the reverse shock, little of it will be directly 
observable initially. The column density of the cool gas also thins out as the supernova ages 
and expands in scale. In such a situation the total luminosity of the reverse shock may 
contribute appreciably, or even dominate, to the bolometric luminosity and is defined as: 

3 



' rad 



1.6 X 10 




erg s 
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Since Kev oc the time dependence of luminosity in cooling case is oc t~^/^^~'^\ 

In Table [3], we tabulate the the timescales up to which the reverse shock remains radiative 
under different conditions. We take the case with density index of 7, 12 and 20 for three 
compositions: solar, helium and oxygen. As one can see from the table, the deeper in the 
ejecta the reverse shock (where the composition is dominated by heavier elements) is, the 
longer time do the reverse shock remain radiative. 

It has been argued that SN 1993J like objects with their high mass loss rate will have ra- 
diative reverse shocks (where cooling is important) even at late epochs (> 100 days) whereas 
SNe with low mass loss rates such as SN 1 999em (a type IIP SN ) with long cooling times will 
develop adiabatic reverse shocks early on (INymark et al.ll2006l ). However in some cases adia- 



batic a nd radiative rey e rse sh ocks co-exist as has been seen in SN 1987A (iGroningsson et al. 



20061 ). iNymark et al.l (120091 ) has discussed that if the ejecta or the GSM is clumpy, the 
adiabatic reverse shock may give rise to slow moving oblique shocks which are radiative in 
nature. 

Soo n after its discovery SN 1 993 J was observed and detected in the X-rays by the 
ROSAT (jZimmermann et al.l 119931 ). Since then, the supernova h as been observed at vari 



ous epochs with a number o f X-ray telescopes, including ROSAT (IZimmermann et allll994 



Immler et al l200lh . ASCA flTanakalll993 



Kohmura et al. 



1994 



Uno et al. 



2002 ) , Cha ndra 



(ISwartz et alll2003l . and this paper), XMM- Newton (jZimmermann fc Aschenbachl l2003l ) and 
Swift (this paper) satellite missions. The Oriented Scintillation Spectrometer Experiment 
(OSSE) on the Compton Observatory had detected the SN in 50-150 keV en ergy band on 



day 1 2 and 30 with very high X-ray luminosity (5 x 10^° erg s ^) on day 12 (ILeising et al. 



19941 ). However, they found that the 50-150 keV emission faded below detection by day 108. 



Suzuki et al.l (119931 ) did a detailed analysis of the early X-ray emission from the SN. They 



carried out hydrodynamical modelling of the collision between the ejecta and a GSM created 
by steady winds {pcsM oc r~^), and claimed that the observed features of X-ray emission can 
be accounted for with thermal free-free emission from the shock heated ejecta. They also 
predicted that the expansion velocities of the ejecta must be high and the density gradient of 
th e ejecta shallow. Later th ey extended their analysis to X-ray observations of first 600 days 



m 



Suzuki fc Nomoto Jl995h . In this paper they assumed a more realistic ejecta model, a 



steep density gradient in the outermost layer followed by a much shallow density gradient in 
the inner H-rich envelope. The density profile suddenly increases and becomes steeper again 
at the interface of H-rich envelope and the He-core. The GSM was assumed to be clumpy 
and cooling effects were taken into account with formation of a dense cooling shell between 
the reverse and the forward shock. Early X-rays from the reverse shock were absorbed in 
this cool shell and X-rays were mainly the thermal emission from the shocked GSM. Their 
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results indicated that the CSM has density distribution of the nature pcsM oc r~^-^ instead 
of PcsM oc r~^, indicating that the CSM was created by non-steady winds. According to 
their model, the CSM has to be more clumpy in the outer layers to explain the observed 
X-ray features. The collisions between the forward shock and the clumps may give rise to 
soft X-rays. Such collisions would take place with many clumps at various distances, and the 
sum of resulting X-ray emissions may form a more or less continuous light curve. In their 
model, all the soft X-rays in first few tens of days were coming from the clumpy CSM, and 
had no contr ibution from the reverse s hell as they were totally absorbed by the cooled shell. 
However, in iNomoto &: Suzukil (119981 ) . they revised their prediction and claimed that the 



cooling dense shell formed between the reverse and the forward shocked shell is clumpy due 
to the Ray leigh- Taylor instability and thus some early X-rays may also leak out from reverse 
shock. They also predicted the possible increase in the soft X-rays around day 500 to be due 
to the density jump from the H-rich envelope to He-layer, modified by the Rayleigh- Taylor 
instabilities. 



Fransson et al.l (119961 ) have studied the early time X-ray emission from SN 1993J in 



detail. Using radio and X-ray observations from the first two weeks, they estimated the mass 
loss rate of ~ 4 x lO~^M0yr~^. Fitting the early radio data using the free-free absorption 
model, they suggested a wind density profile: pcsM oc r~^-^. They did the detailed analysis 
of X-ray data from ROSAT, ASCA and OSSE on the Compton Observatory. According to 
their conclusion, the early hard X-rays arise mainly from the CS shock, whereas the high 
column density of the cooling shell between the CS shock and the radiative reverse shock 
absorbs most of the soft X-ray emission arising from the radiative reverse shock. After about 
50-100 days the column density of the cooling shell decreases enough for the X-rays from 
the reverse shock to leak out and dominate the spectrum. They also found that the X-ray 
observations are best fit with the flatter CSM density profile of pcsM r"^-''. However, in 
their later work iFransson fc Bjornssonl (119981 ) claimed that the CSM density profile is more 
sensitive to the radio observations and they found the radio observations to be fitting well in 
PcsM OC r~^ profile. The earlier discrepancy in radio band was claimed to be due to neglecting 
the synchrotron-self ab sorption process which play s an important role in the radio absorption 
for SN 1993J. However, iMioduszewski et all (120011 ) included both synchrotron self-absorption 
and free-free absorption in their analysis and claimed that the radio emission could not be 
fit with a CSM declining as r~^, but that a density profile going as r~^'^ provided a much 
better fit. In this paper, we use the CSM profile to be pcsM oc r~^, though we explore the 
possibility of different profiles in §61 
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4. Observations and Analysis 

As mentioned in ^ SN 1993J lias been observed at various epochs with multiple X-ray 
telescopes since its discovery. We describe below these observations in detail. We analyzed 
the 2005 and 2008 Chandra data and the Swift data at several epochs between 2006-2008. 
We also reanalyzed the 2001 XMM-Newton and 2000 Chandra datasets and extracted un- 
absorbed fluxes in 0.3-2.4 keV, 2-8 keV and 0.3-8 keV bands. For ROSAT and ASCA, 
we extracted unabsorbed fluxes in the above bands using best flt parameters mentioned in 
various references (see below). In our flts, the column density is a measure of the absorbing 
cool shell in the SN plus the galactic absorption shell. The column density was set to zero to 
extract the absorption corrected luminosities. Below we describe the various SN 1993J X-ray 
observations in detail. The details of all observations used in this paper are are summarized 
in Table [TJ 



4.1. ROSAT observations 



X-ray observations of SN 19 93J by i? 05"^ T between days 6 to 1181 after explosion was 
reported by llmmler et al I (120011 ). in which both PSPC and HRI observations were included. 
HRl data were binned in observation blocks of length 5-20 ks exposure time whereas the 
PSPC observation blocks had integration times ranging 2-18 ks long. They had assumed a 



consta nt absorption column density and extracted the fluxes. IZimmermann &: Aschenbach 



( 120031 ) reanalyzed all the ROSAT data and reported the observations up to day 1800. They 
flt all the spectra assuming column density as a free parameter and with a thermal plasma 
model with a single temperature (vmekal in XSPEC) flxing the element abundances to 
that obtained from, e.g. XMM-Newton PN spectrum. They found that due to the high 
temperatures in the early ROSAT observations there is almost no difference between fltting 
either with solar or with the XMM-Newton elemental abundances. We use the best flt 



parameters quoted in Table 3 of IZimmermann fc AschenbachI (120031 ) to extract the ROSAT 
fluxes of the SN at various epochs. Table [2] shows the unabsorbed luminosities for the 
ROSAT-HRI and ROSAT-PSPC observations in the energy range of 0.3-2.4 keV. 



4.2. ASCA observations 



ASCA flux po ints were extracted from iKohmura et al.l (119941 ) , lUno et al.l (12002! ) and 
Swartz et all (120031 ). The large Field of View (FoV) of ASCA showed the presence of the SN 
host galaxy and a bright X-ray binary source at respectively 3' and 1' away from the SN. 
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Since ASCA PSF has resolution of 3', t he bright X-ray b inary was s trongly contamiii ating; the 
SN flux. To avoid this contamination iKohmura et al.l (119941 ) and lUno et al. isolated 
SN 1993J flux through a one- dimensional image fltting, by using 1' x 2' rectangular boxes 
containing the SN and the X-ray binary separately, and by modelling the one- dimensional 
intensity proflle. They claimed to get rid of almost all the contamination from the M81 and 
the X-ray binary source, by this method. 

To estimate more reliable ASCA fluxes of SN 1993J, we attempted to utilize the late 
Chandra observations. Due to its high angular resolution, Chandra is able to easily separate 
the SN, and the nearby X-ray binary and M81 nucleus. Our aim was to extract the uncon- 
taminated flux of the X-ray binary and then subtract it from all the ASCA observations. 
However, for this technique to work, we needed to make sure that the X-ray binary is not 
a time variable source. To determine this, we analyzed two Chandra datasets, one of 2000 
May 7 and another on 2005 Jun 01 and extracted the flux of the X-ray binary from both 
the data respectively. Our analysis shows that the X-ray binary is a highly variable source, 
especially in the hard X-rays. The 0.3 — 8.0 keV flux decreases by a factor of two in flve 



years, from 2.79 x 10~^^ erg cm~^ s~^ in May 2000 to 1.40 x 10~^^ erg cm~^ s~^ in June 
2005. The change in the SN flux at 2 — 8 keV is by more than a factor of three in the two 
observations (2.26 x 10~^^ erg cm~^ s~^ in May 2000 to 6.9 x 10~^^ erg cm~^ s~^ in June 
2005). Thus, in view of high variability of the X-ray binary source, it is not advisable to 
reanalyze all the ASCA datasets to remove the contribution of the X- ray binary on the basis 
of Chandra flux of th e binary. We use the spectral flt parameters of iKohmura et al.l (119941 ) 
and lUno et al.l (120021 ) to convert the count rates into unabsorbed fluxes in bands 0.3-2.4, 
2-8 and 0.3-8 keV fluxes (Table [2]). 



4.3. Chandra observations 



Chandra flrst observed SN 1993J on 2000 March 21 and 2000 May 07 under ObsIDs 
390 and 735 with ACIS-S bv ISwartz et~a3 J2003h . They showed that SN 1993 J faded since 
its discovery and its spectrum softened. At this time SN 1993J was ~ 2600 days old and 
displayed a complex thermal spectrum from a reverse shock rich in Fe L and highly ionized 
Mg, Si, and S but lacking Oxyge n. The unabsorbed luminosity reported in Table [2] are 
extracted from lSwartz et all (120031 ). 



Chandra ACIS-S observed the SN starting from May 26, 2005 to July 6, 2005 on 15 
different occasions (Obs Ids: 5935-5949; PI: Pooley). Each observation was around 11-12 
ks. We analyzed this data using CIAO analysis threads and XSPEC. Event 2 flies (pipeline 
processed flies) were used for the data analysis. Standard methods were used to analyze the 
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data (IChandra et al.ll2005l ). We combined the whole data set in two groups to increase the 
signal to noise ratio and extracted the unabsorbed fluxes at various epochs in 0.3 — 2.4 and 
2.0 — 8.0 keV ranges. We also did spectral analysis of ~ ISOks Chandra ACIS-S data in 2005 
given its large exposure. We get an excellent fit to the data using a two- component thermal 
plasma model (reduced of 1.09) with absorption column density of (6.0 ±1.5) x 10^° cm~^. 
This is close to the Galactic absorption in the direction of M81, which is expected at such late 
epochs when the cool shell absorption may have become insignificant. The two temperatures 
obtained from the best fit are 0.73 ± 0.04 keV and 2.21 ± 0.24 keV. Other models do not fit 
the data well. The Non-equilibrium ionization (NEI) model gives reduced = 1-88 for 104 
degrees of freedom and power- law model gives reduced = 3.48 for 105 degrees of freedom. 
The two temperature bremsstrahlung model yields a reduced = 1-90 for 103 degrees of 
freedom. Hence, the two-component thermal plasma model is the most plausible one. We 
could not fit the various lines to the data given its sparseness. We plot the spectral fit to 
this data in Fig. [H 

We also extracted the archival data observed with Chandra ACIS-S in HETG mode at 
15 occasions from Feb 24, 2005 to Aug 12, 2006 (Obs IDs: 6174, 6346-47, 5600-01, 6892- 
6901; PI: Canizares) and analyzed it. These observations were taken for M81, the host 
galaxy of SN 1993J which is 3' away from the SN. SN 1993J fiux was a byproduct of these 
observations. Half of the observations were centered on ACIS chip S2 and the rest of the half 
on S3. We combined all the datasets in these two major sets and extracted the fiux using 
best fit models. In this case, we fixed the column density to be 0.06 x lO^^cm"^, the one 
obtained from the 2005 Chandra observations (see above). 

Latest observations of SN 1993 J with Chandra was taken on 2008 Feb 01 (PI: Immler). 
The observations were taken in ACIS-S mode without any grating. The total exposure time 
was 14.80 ks, out of which we could extract 10 ks of good data. The count rate in this 
observation was (2.22 ± 0.15) x 10"^ cts in full ACIS band. We converted these count rate 
into fiux using Nh = 0.06 x lO^^cm"^ and temperature of 0.7 keV. 



4.4. XMM-Newton observations 



XMM-Newton observed SN 1993J on 22nd April 2001 around eight years after the 



explo sion for about 132 ks duration, under obsID 0111800301 (jZimmermann fc Aschenbach 



20031 ). Data from the EPIC P N camera, run in small window mode , and the M0S2 camera in 
imaging mode were obtained. IZimmermann &: AschenbachI (120031 ) did the detailed analysis 
of the data and fit the X-ray spectrum with 2-component thermal plasma model. They 
claimed emissions from highly ionized Mg, Si, S, Ar, Ca and complex Fe. We used their best 



- 10 - 



fit parameters to extract the unabsorbed luminosity in 0.3-2.4 keV, 2-8 keV and 0.3-8 keV 
bands. 



4.5. Swift observations 



The X-Ray Telescope (XRT lBurrows et a on board the Swift Observatory (IGehrels et all 



20041) observed SN 1993J on 2005 April 21, 2005 August 25, 2006 June 24, 2006 November 



18 — 20 and on 2008 January 09-15. We combine the data closely spaced in time, i.e. 2005 
April 21-August 25, 2006 June 24-November 18, and 2008 January 09-15, to increase the 
sensitivity of the data. 

The HEASOFtQ (v6.2) and Swift software (v2.6.1, build 20) tools and latest calibration 
products were used to reanalyze the data. X-ray counts were extracted from a circular region 
with a 10 pixel (24") radius centered on the optical position of the SN. The background was 
extracted locally from a source-free region of 30" radius and corrected for the 100% encircled 
energy radius, to account for the detector and the sky background, and for residual diffuse 
emission from the host galaxy. SN 1993J is detected in X-rays in the 21 ks XRT observation 
from 2006 November 18 — 20 at a 7.0a significance of source detection and a net count rate 
of (4.1 ± 0.6) X IQ-^ cts s"^ (0.2-10 keV). Table [2] gives details of these observations. 



5. Results and Interpretation 

We tabulate the details of all the above observations described in fg] in Table [TJ We 
report the unabsorbed luminosities in 0.3-2.4, 2-8 and 0.3-8 keV bands in Table [21 We plot 
the light curves constructed in these bands from the above observations in Fig. [2l The light 
curve in 0.3-8 keV shows linear decline with powerlaw index of —0.65 for first few hundred 
days. At late epochs, from day 2500 onwards, the light curve declines at t^^. Due to the lack 
of observations between day 500 to 2500, it is not possible to predict its evolution in this 
time range. However, the fiux seems to have increased in this gap of ~ 2000 days, contrary 
to the expected decline. The light curve in 2-8 keV seems to have declined much faster 
with a powerlaw index of —1 throughout the all the observations. In this band too, there 
seems to be a possible increase in the luminosity in the gap of 500-2500 days. Fortunately, 
ROSAT-PSPC observations have covered this gap with many data points in the softer band. 
The light curve in 0.3-2.4 keV shows an overall decline in the light curve with a powerlaw 



http: / /heasarc. gsfc.nasa.gov/docs/software/lheasoft/ 
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index of —0.25 for until about 1500 days. Thereafter, the hght curve is consistent with 
dechne. However, after day ~ 200 — 300 the hght curve reveals a sudden drop in the 
luminosity which rises back slowly and reaches back to the —0.25 index profile by day 800. 
There is another possible small dip around day 1000-1500. 

As discussed in ^ the X-rays may originate either from the forward shock or from the 
reverse shock. Since the density and temperature conditions of the two shocks are widely 
different, this information can help us identify the X-ray emission regions in a given energy 
band. The CS temperature can be written as 

r,, = 1.36xlO»f^Vf-:.J^VK. 



n-2j VlO^kms 

Radio VLBI of S N 1993J suggests an expansion velocit y of the order of 7000 — 8000 km s~^ 



at current epoch (IBartel et all 120071 : iBartel et al.l 120021). H-a observa tions suggest that the 



ejecta velocity has to be at least 7500 km s~^. iFransson et al.l (Il996l ) have argued that the 
value of density power law index n is 8 — 12. Putting these values in the above equation 
implies that CS temperature should be at least Tcs > 50 keV. At such high temperature of 
the CS shock, there will be very little emission which will come out in the 0.3-8 keV band. 
At early epochs the velocities would be much higher resulting in high Tqs, of the order of at 



l east f ew tens of keV and thus will have little flux in the energies below 8 keV. iFransson et al. 



(119961 ) have also indicated that all the flux below 10 keV is mostly arising from the reverse 
shock on the basis of early X-ray observations of the SN. The above discussion indicates 
that the 0.3-2.4 keV and 2-8 keV bands are at much lower temperatures to account for the 
emission from the forward shock; and most probably the X-rays in these bands are arising 
mostly from the reverse shock. Our claim is further strengthened by the hydrodynamic 
simulation (see $ED, which also seems to suggest the reverse shock origin. 

X-ray light curve behaviour can give information about the nature of the reverse shock 
from which majority of X-rays is coming. The adiabatic shock luminosity declines with a 
powerlaw of index of —1, whereas the radiative reverse shock declines much slowly. The slow 
decline of 0.3-8 keV light curve at early epoch in consistent with the radiative nature of 
the reverse shock at these epochs. The late time —1 decline is consistent with an adiabatic 
reverse shock. Due to the lack of data it is not possible to pin point the transition from 
radiative shock to adiabatic shock. However the soft band light curve without this gap of 
2000 days clearly indicates that the shock made the transition from a radiative reverse shock 
to an adiabatic on around day 1500-2000. Initially hard band emission is much higher than 
the soft band emission but the latter dominates from day 200 onwards. Since SN 1993J 
emission is becoming softer with time, this kind of luminosity evolution is expected. The 
fast decline of luminosity in 2-8 keV band also indicates that very little of this radiation is 
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being absorbed by the cool shell in this band and most of it is coming out as adiabatic shock. 
The overall light curve evolution is consistent with an early radiative reverse shock turning 
into an adiabatic s hock after around 1000 days. This kind of behaviour has been predicted 



by iNymark et al.l (120061 ) for such supernovae. The sudden decline and gradual rise of the 
X-ray flux between day 200-800 and a possible dip around day 1200-1500 may be due to 
the density fluctuations in the ejecta. ISuzuki fc Nomotd (Il995l ) has predicted a possibihty 
of a jump in the SN 1993J light curve. They claim that while moving backwards into the 
ejecta, the reverse shock encounters a density jump at the interface of the hydrogen envelope 
and the Helium core. This density jump at the interface of th e two regions will show up as 
luminosity jump in the light curve. iMioduszewski et all (120011) clearly show the effect of this 
density jump on the radio emission. ISuzuki &: Nomotd (119951) have calculate d the p ossible 



time for this jump to be around 500-1000 days, although IMioduszewski et all (120011 ) find it 
to be slightly la ter. It is probable that w e are seeing this effect in the light curves of SN 
1993 J. However, ISuzuki fc Nomotd (119951 ) model is 1-D and the strong density gradients are 
likely to be smoothed by instabilities. 

We use the cooling time expression in Eq. |2] to derive a timescale for which a reverse 
shock may remain radiative for the relevant parameters for SN 1993J mentioned in previous 
sections. A rapid cooling of the gas ahead of the shock is required such that at any stage 
the cooling timescale is less than the elapsed timescale so that a layer of cool and dense 
absorbing gas forms, i.e. tcooi/^ < 1- This condition along with Eq. |2] gives 



t < 



(n-3)(n-4)(n-2) 
605 



3.34 




10^ kms-1 



-5.34 



days. 



(4) 



For n = 12, this equation gives t < 260(M_5/Mwi)(Vej/10^kms"-^)"^-^^ days. We tabulate 
radiative timescales for various values of Vej and M_5 in Table |3l This indicates that it is 
possible for reverse shock to remain radiative at late epochs, as indicated by our light curves. 



6. Hydrodynamic simulation and computed X-ray light curves 

The analytic calculations in this paper assume a CS density profile that varies as r~^. 
However, in numerical simulations, we do not constrain ourselves with and explore a 
larger parameter space. We have carried out several analytic simulations to find the best 
possible CSM density structure that provides an adequate fit to the observed X-ray data. We 
chose to reproduce the data in the hard X-ray band, because this band is likely to behave 
adiabatically (as evident from our 2-8 keV light curve) and hence h ave no complicated 
line emission due to the cooling effects. Our simulations are based on IMioduszewski et al 
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( 120011 ) for a spherica lly symmetric system usi ng high resolution VH-1 three dimensional 
finite difference code. iMioduszewski et all (1200 ll ) produced simulated radio light curves using 
a detailed radiative transfer calculation and fit the radio light curve of SN 1993J. These 
hydrodynamic results were updated in iBartel et all (120071 ) to match updated VLBI data 
from 1993J. Herein we use an updated version of this hydrodynamic calculation to simulate 
the observed X-ray light curves of SN 1993J reported in this paper. The code computes the 
interaction of the ejecta with a CSM whose density profile leads to the calculation of X-ray 
fiux density of the SN at each epoch. The calculations were carried out in one dimensional 
Lagrangian coordinates and remapped on to an Eulerian grid. The shocked interaction 
region between the forward and reverse shocks was adequately resolved, which is important, 
as most of the X-rays are emitted in this region. 

The hydrodynamic runs use the 4H47 ejecta density distribution model of Shigeyama 
and Nomoto (courtesy of K. Nomoto 1999, private communication). This model had an 
ejecta mass of 3.12 Mq, the pre-SN progenitor star radius of 350 Rq and an envelope mass 
of O.47M0. The mass fraction of helium in the envelope was about 0.79. The explosion 
kinetic energy was 10^^ erg. The surrounding density profile has a slope that varies over 
time, from r~^'^ in the very early stages, to r~^-^ up to 10^'' cm, and then a steeper drop 
to r~^-^ at late epochs. This CSM density profile was used to adequately fit the VLBI 
data in iBartel et all (120071 ) . and produced radii and velocities comparable to the radio and 
optical (Ha) observations. The computation of radiative emission from the hydrodynamic 
model further requires the electron temperature, whereas the the hydrodynamics gives the 
post shock temperature, which is the temperature of the ions Tion- We have found that 
an electron temperature about 15 % of the ion temperature provides an adequate fit to 
the emission. This is a somewhat simplistic assumption. We calculated the X-ray emission 
between 1-6 A (2-10 keV). The observational analysis show that the absorption by the cool 
shell is less significant for the hard X-ray emission, thus we ignore it. Line emission may 
also be ignored in this hard band. We used abunda nces that ar e appr o priate for the a r nbient 
medium and SN ejecta in SN 1987A, as outlined by iLundqvistI (119991 ). iNymark et al.l (120091 ) 
also found the SN 1987A abundances to be a reasonable set of abundances for the SN 1993J. 



The io nization of various species were computed by a model due to IShuU fc van Steenberg 
(119821 ). Differences between various ionization models are small and do not alter the hard 



X-ray fiux significantly, or the fundamental result that the X-ray fiux arises from reverse 
shocked ejecta. 

Given the density and temperature in every zone, the hydrodynamic mode l calculates 



the fr ee-free and bound-free luminosity from that zone using the CHIANTI code ( jPere et al. 



19971 ). This is being done for 10000 zones, over each time-step. The contribution from each 
zone is added up to give the total luminosity. Having the contribution from each zone also 
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allows us to determine exactly the part of the density profile that is contributing most to 
the X-ray emission. The computed X-ray luminosity in the hard band is shown in Fig. [31 
Looking at the simplicity of our run, the match between the simulated light curve and actual 
luminosity is encouraging. Between 10 and about 400 days the luminosity decreases only 
slowly. Some of the later peaks are due to the reverse shock running into ejecta structures 
which compresses them. After 1000 days, the slope becomes steeper (L oc t~^), in accordance 
with the data. We find that the hard X-rays arise from just behind the reverse shock. The 
region of hard X-ray production expands inward in the ejecta over the next decade or so 
till they begin to arise from most of the shocked ejecta, between the reverse shock and 
the contact discontinuity. After about 50 years, the forward shock starts to dominate the 
emission. This confirms our prediction that almost all the X-rays below 8 keV over the 
observed time period are coming from the reverse shock. 

Our hydrodynamic simulation models are carried out in spherical symmetry. It is en- 
couraging that even using these simple hydrodynamic models it is possible to understand 
qualitatively the evolution of the SN, and explore several basic features, such as the fact that 
the hard X-ray emission is coming predominantly from the reverse shock. 



7. Discussion 



7.1. Column Depth and the cool shell 



Immler et al I (120011 ) had analyzed ROSAT data fixing the column depth to be 4 x 



10 cm , same as the Galactic absorption c o lumn density. However, the reanalysis of 



the same data by IZimmermann fc Aschenbachl (120031 ) and letting column density to be a 
free par ameter rev e aled that column densities are much higher than the Galactic column 
density. lUno et al.l (120021 ) fit the ASCA data with two component thermal plasma model and 
found that the low temperature component has much higher column density than the high 
temperature component. These high column densities in excess to the Galactic absorption 
can be attributed to the absorption by an additional cool shell. Since most of the X-ray 
emission below 8 keV is coming from the reverse shock, this indicates that both the ROSAT 
and the ASCA data have revealed the presence of a cool shell between the forward and a 
reverse shock, thus presence of a radiative reverse shock. However, the column density seems 
to reduce with time and after around 2000 days, the best fit column densities are consistent 
with that of the Galactic absorption. This indicates that the reverse shock has most likely 
become adiabatic by this epoch. 
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7.2. Hardness Ratio and electron temperature 

In Fig. m we plot the ratio of SN 1993J luminosity in the 2-8 keV versus 0.3-2.4 keV 
bands. The figure clearly shows that the hard X-ray emission dominates for first ~ 100 — 200 
days. After ~ 200 days, soft X-ray emission starts to take over and continues to dominate. 
At current epoch, i.e. 15 years after the explosion, the soft X-ray is dominant by around an 
order of magnitude. 

Since most of the fiux below 8 keV is coming from the reverse shock, the ratio of the 
reverse shock free- fre e luminosity in 2-8 keV band to 0.3-2.4 keV band can be written as 



(jFransson et al.lll996l ): 



^ ^ -^2.0-8.0keV ^ /2^0keV ( ^ ) "'^ exp ( - ) d-E 

Here {E/kT)~^'^ is the Gaunt factor in the given energy regime. This equation puts some 
very interesting constraint on the temperature of the reverse shock. The current hardness 
ratio of 0.14 corresponds to the electron temperature of 1.05 keV. This equation also shows 
that electron temperature was ~ 1.5 keV around day 500 and it has evolved very slowly 
since then to reach the current value of 1 keV. Since free-free emission dominates above 



temperature 2 keV (iFransson et al.l Il996l ). our derived temperature of 1 keV at current 
epoch will contribute mostly to the free-free continuum emission, as well as possibly to the 
line emission. 



For maximum electron temperature of ~ 30 keV (iFransson et al.lll996l ). the hardness 
ratio is 1.45. In fact, the Eq. [5] above reaches asymptotic upper limit of ~ 1.5 which is 
around day 100, irrespective of temperature. However, this demonstrates that in the initial 
X-ray observations, where hardness ratio is larger than 1.5, the contribution between 2 — 8 
keV is not solely because of the reverse shock but may have some contribution from the CS 
shock too. After day ~ 100, all the X-ray fiux below 8 keV comes from the reverse shock 
(see Figure H]). 



7.3. Ha and the X-ray luminosity evolution 

Ha emission in SNe arises initially by radioactivity and later by reprocessing of X-rays 
produced due to ejecta-wind interaction. In case of adiabatic shocks, the Ha emission may 
come from the unshocked material heated by the X-rays coming from the adiabatic shock 
waves. If the reverse shock is radiative, then the dense cool shell between the reverse and 
the forward shock may also give rise to significant Ha emission. 
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There have been extensive Ha observations of SN 1993 J. iMatheson et al.l (2000a]Jbl ') 



Matheson et al 



J2000bf i 



describe the detailed optical spectra of SN 1993 J up to 2500 days. In 
they also measure the l inewidths o f boxy Ha profile, an indication of the CS interaction, 
from day 433 on-w ards. iPatat et all (119951 ) reports Ha fluxes between day 171 to day 367. 
Zhang et al.l (12004 ') repor t photometric observations of SN 1993J from the year 1995 to 2003. 
Houck fc FranssonI (119961 ) showed that the decay of ^^Ni provides the ultimate power source 
for Ha emission up to 250 days. After 250 days, continued presence of strong broad lines of 
Ha in SN 1993J have been explained as a result of photoionization by X -rays and U V emissio n 
from the radiative reverse shock propagating into the supernova ejecta (IFransson et al.ll2005l ). 
Pat at et all (119951 ) show that late time Ha emission from SN 1993J can be described by a shell 
of hydrogen between 7500 to 11,400 km s~^. Similar velocity range for emi tting gas shells 



for b road, box-shaped UV lines are seen in the HST spectra of SN 1993J (jPransson et al. 
20051 ) apparently coming from an ejecta and a cool dense shell. 



We tabulate these Ha luminosities after day 170 in Table HI We also plot the Ha 
luminosities with soft and hard band X-ray luminosities in figure [51 The figure shows that 
the Ha luminosity seems to follow the hard X-ray band luminosity evolution after a few 
hundred days, but with significant efficiency. However, it fails to trace the soft band X-rays 
completely. The significant rate of conversion of the kinetic luminosity of the ejecta wind 
interaction into broad Ha emi ssion by the repro cessing of the X-ray lumin osity of the r everse 
shock wave was also noted by lPatat et all (119951 ). It should be noted that lPatat et all (Il995l ) 
have cautioned that possibly up to 30% of the emission at Ha wavelengths in SN 1993J may 
originate from an unidentified band of emission around 6600 A which is hypothesized due 
to a broad blend of emission lines extending between 6050-6800 A as seen to be present in 
type Ib/Ic SNe. 

The Ha seems to trace the hard band evolution but the efficiencies are 30-50% of the 
X-ray production in this band, whereas the expected efficiency is typically no more than 
1-5% of the X-ray luminosity. This indicates that there could be another component which 
is contributing significantly to Ha fiux. The clumps in the ejecta can be one such candidate 
which can give rise to Ha with high efficiency. The non-smooth Ha evolu tion also probably 
indicates towards the possibility of significant Ha origin from the clumps. iPatat et all (119951 ) 
also invoked an earlier suggestion by lChugail (119931 ) where the dumpiness of the wind material 
and possibly the clumpy structure of the dense and Rayleigh Taylor unstable region in 
the ejecta helped a t tain a more efficient transformation of kinetic energy into radiation. 
Chugai fc Danzigeil (119941 ) proposed a model for SN 1988Z, in which a radiative shock wave 
is driven into the dense cloud by thermal and dynamical pressures behind the main (blast 
wave shock and reverse shock) shock waves or by the dynamical pressures of the expanding 
unshocked ejecta. The shock in the cloud is significantly slower due to the higher density 
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in the cloud and the shocked gas cools by soft X-ray or UV emission which then pumps the 
optical emission from the cool dense material behind the radiative shock wave in the cloud. 



8. Conclusions 

In this paper, we have presented the complete X-ray light curves of SN 1993J in 0.3- 
2.4 keV, 2-8 keV, and 0.3-8 keV bands. We demonstrate that most of the emission below 
8 keV comes from the reverse shock, except for very early soft band emission when the cool 
shell absorbed all the soft flux from the reverse shock. The light curves reveal that the 
reverse shock is radiative for around initial 1000 days and then it becomes adiabatic at later 
epochs, as expected for such supernovae jNvmark et al.ll2006h . The evolution of column 



density too seem to indicate this. The column density was higher than the Galactic column 
density during the ROSAT, ASCA and early Chandra and XMM observations, and later 
became comparable to the Galactic absorption, indicating presence of a cool shell during 
around first 1000 days. We demonstrate that for SN 1993J parameters, it is possible for the 
shock to be radiative at such late epochs. 

We have carried out numerical hydrodynamic computations, and calculated the hard X- 
ray flux from the same, which agree reasonably well with the observed data. Our simulations 
clearly show that all the emission below 8 keV is indeed coming from the reverse shock, 
strengthening our claim. 

The large fraction of Ha flux in comparison to the X-ray emission indicates possibility 
of clumps in the ejecta. 

One of the important questions is how these young supernovae evolve towards a Super- 
nova Remnant. Since young SNe and older SNRs are both products of explosive events a 
natural question is whether the former evolve continuously into the latter with time, whether 
similar emission mechanisms differing only in length-scales and time-scales operate or if the 
SNe fade away only to switch on later as SNRs with a different mechanism of radiation. One 
of the very well studied SNR in our galaxy, Cassiopeia A was thought to be of the Type lib 



or Iln IChevalier fc Oishil (120031 ). However, it h as recently been co nclusively classified as a 



Type lib SN, the same type as that of SN 1993 J (IKrause et al.ll2008l ). Late time observations 
of SN 1993J may provide an important link between young SNe and SNRs. SN 1993J had a 
large X-ray flux, which, coupled with it being a nearby SN, made it easily observable for a 
long time with a well-sampled light curve. This has made SN 1993J one of the best studied 
extragalactic supernovae, and provided a wealth of information on SN evolution and their 
circumstellar interaction. Further observations of this SN for a long time is sure to provide 
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much valuable information on these questions. 

We thank the anonymous referee for useful comments. P.C. thanks Roger Chevalier for 
useful discussions. We thank Richard McCray for discussions and his comments on an earher 
draft of this manuscript and thank Claes Pransson, Roger Blandford, Alex Filippenko, Subir 
Sarkar, Stephen Smartt for several illuminating discussions. We embarked on this work at 
the 2007 Aspen Workshop on "Supernova 1987A: 20 Years After" and wc thank the organiz- 
ers of that Workshop. A.R. thanks the participants of the Darjeeling School and Workshop 
on Supernovae and GRBs in May 2008 where this work was discussed and Manjari Bagchi 
and Sayan Chakraborti for their technical assistance. P.C. is a Jansky fellow at National 
Radio Astronomy Observatory. The National Radio Astronomy Observatory is a facility 
of the National Science Poundation operated under cooperative agreement by Associated 
Universities, Inc. At Tata Institute this research was supported by the Eleventh Pive Year 
Plan Project No. llP-409. V.V.D. received support from award 7^ AST-0319261 from the 
NSP. This research has made use of data obtained through the High Energy Astrophysics 
Science Archive Research Center On-line Service, provided by the NASA/Goddard Space 
Plight Center. We thank the Chandra X-Ray Observatory team for carrying out the obser- 
vations for Peb 2008 under guest observing program. CHIANTI is a collaborative project 
involving the NRL (USA), RAL (UK), MSSL (UK), the Universities of Florence (Italy) and 
Cambridge (UK), and George Mason University (USA). 

REFERENCES 

Bartel, N., Bietenholz, M. P., Rupen, M. P., Dwarkadas, V. V. 2007 ApJ668, 924 
Bartel, N., et al. 2002, ApJ, 581, 404 

Burrows, D. N., et al. 2005, Space Science Reviews 120, 165 

Chandra, P., Ray, A., Schlegel, P.M., Sutaria, P. K., & Pietsch, W. 2005 ApJ, 629, 933 
Chevaher, R., & Oishi, J. 2003 ApJ, 593, L23 

Chevalier, R., & Pransson, C, 2003, Supernovae and Gamma-Ray Bursters, (Ed. K. Weiler), 
Lecture Notes in Physics, 598, 171 

Chevaher, R. A., & Pransson, C. 1994, ApJ, 420, 268 

Chevaher, R., & Soker, N. 1989 ApJ341, 867 



-19- 



Chevalier, R. 1982 ApJ, 258, 790 

Chugai, N. 1993 ApJ, 414, LlOl 

Chugai, N. & Danziger, J. 1994 MNRAS, 268, 173 

Dere, K.R et al, A&A Suppl. Ser, 125, 149 

Filippenko, A. V., & Matheson, T. 2005, lAU CoUoq. 192: Cosmic Explosions, On the 10th 
Anniversary of SN1993J, 37 

Fihppenko, A. V., Matheson, T., & Woosley, S. E. 1993, lAU Circ, 5787, 1 

Fransson, C, et al. 2005, ApJ, 622, 991 

Fransson, C, & Bjornsson, C.-I. 1998 ApJ, 509, 861 

Fransson, C, Lundqvist, R, & Chevalier, R.A. 1996 ApJ, 461, 993 

Fransson, C, & Sonneborn, G. 1994 Frontiers of Space and Ground Based Astronomy, ed. 
W. Wamsteker et al. Astrophysics an Space Science Library, 187, p. 249 

Freedman, W. L., et al. 1994, ApJ, 427, 628 

Gehrels, N., et al. 2004, ApJ 611, 1005 

Groningsson, R, Fransson, C, Lundqvist, R, Nymark, T., Lundqvist, N., Chevalier, R., 
Leibundgut, B., & Spyromilio, J. 2006, A&A, 456, 581 

Houck, J. C. & Fransson, C. 1996 apj 456, 811 

Immler, S., Aschenbach, B. & Wang, Q. D., 2001, ApJ, 56, L107 

Krause, O., Birkmann, S. M., Usuda, T., Hattori, T., Goto, M., Rieke, G. H., & Misselt, 
K. A. 2008, Science, 320, 1195 

Kohmura, Y., et al. 1994, PAS J, 46, L157 

Leising, M. D., et al. 1994, ApJ, 431, L95 

Lundqvist, R 1999, ApJ, 511, 389 

Matheson, T., et al. 2000, AJ, 120, 1487 

Matheson, T., et al. 2000, AJ, 120, 1499 



- 20 - 



Matzner, CD., & McKee, C.F., 1999 ApJ, 510, 379 

Maund, J., Smartt, S.J. et al, 2004, Nature, 427, 129 

Mioduszewski, A. J., Dwarkadas, V. V., & Ball, L., 2001, ApJ, 562, 869 

Nadyozhin, D.K., 1985, Ap&SS, 112, 225 

Nomoto, K. et al., 1993, Nature, 364, 507 

Nomoto, K., & Suzuki, T. 1998, The Hot Universe, 188, 27 

Nymark, T. K., Chandra, P. & Fransson, C, 2009, To appear in AfcA. [arXiv:0812.2252l 

Nymark, T. K., Fransson, C. & Kozma, C, 2006, A&A, 449, 171 

Patat, F., Chugai, N., & Mazzah, P.A. 1995, Astron. Astrophys., 299, 715 

Podsiadlowski, P., Hsu, J. J. L., Joss, P. C, Ross, R. R., Nature, 364, 509 

Ray, A., Singh, K. P., & Sutaria, F. K., 1993, J. Astrophys. Astron., 14, 53 

Ripero, J. & Garcia, F., 1993, lAU Circ, 5731, 1 

Shull, J. M., & van Steenberg, M., 1982 ApJS, 48, 95 

Suzuki, T., & Nomoto, K. 1995, ApJ, 455, 658 

Suzuki, T., Kumagai, S., Shigeyama, T., Nomoto, K., Yamaoka, H., & Saio, H. 1993, ApJ, 
419, L73 

Swartz, D. A. et al. 2003, ApJS, 144, 213 

Swartz, D. A., Clocchiatti, A., Benjamin, R., Lester, D. F., & Wheeler, J. C. 1993, Nature, 
365, 232 

Tanaka, Y. 1993, lAU Circ, 5753, 1 

Truelove, J. K. & McKee, C. F. 1999, ApJS, 120, 299 

Uno, S., et al. 2002, ApJ, 565, 419 

Van Dyk, S. D., Weiler, K., Sramek, R., et al. 1994 ApJ, 432, L115 

Woosley, S. E., Eastman, R. G., Weaver, T. A., Pinto, P. A., 1994, ApJ, 429, 300 

Zhang, T., Wang, X., Zhou, X., Li, W., Ma, J., Jiang, Z., & Li, Z. 2004, AJ, 128, 1857 



Zimmermann, H.-U. & Aschenbach, B. 2003, A&A, 406, 969 
Zimmermann, H.-U., Lewin, W., Predehl, P. et al. 1994, Nature, 367, 621 
Zimmermann, H. U., et al. 1993, lAU Circ, 5748, 1 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-22- 



Table 1. Details of X-ray observations of SN 1993 J 
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wCl Zl.lo — WCl LL.OO 


Ac;p A 




■lOOOQOOO 
ozuuyuuu 


40.4 


1995 


Apr 13.71-May 04.90 


ROSAT 


HRI 


RH600740N00 




1995 


Oct 19.00 


ROSAT 


HRI 


RH600881N00 




1996 


Apr 15.88-May 07.83 


ROSAT 


HRI 


RH600882N00 




1996 


Nov 04.50 


ROSAT 


HRI 


RH600882A01 




1997 


Mar 31.10 


ROSAT 


HRI 


RH601001N00 




1997 


Sep 30.44-Oct 16.00 


ROSAT 


HRI 


RH601002N00 




1998 


Mar 26.10 


ROSAT 


HRI 


RH601095N00 




2000 


Mar 21.03 


Chandra 


ACIS-S (NONE) 


390 


2.4 


2000 


May 07.23 


Chandra 


ACIS-S (NONE) 


735 


50.6 
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Table 1 — Continued 



Date of 


Mission 


Instrument 


Observation 


Exposure 


observation 








ID 


(ks) 


2001 


Apr 22.32 


Y MM- Np^H-nn 


EPIC-PN 


01 1 1 800301 


1 39 


2005 


Aug 14.41 


f .H nin n vn 


ACIS-S 


(HETG) 


uuuu 




2005 


Jul 19.60 


( . n nn /I T'/i 

y^liiJbllUjI tb 


ACIS-S 


(HETG) 


uuux 


RA 1 


2005 


Jul 14.81 


(in (IT) nvn 


ACIS-S 


(HETG) 


6S47 


64 7 


2005 


Jul 14.07 


f .11 nti n T/~t 

KJI bill vUjI Uj 


ACIS-S 


(HETG) 


6'^46 


55 2 


2005 


Feb 24.29 


( ,n HTi n TTI 
K^l LUjI LUjI Ui 


ACIS-S 


(HETG) 


61 74 


45 9 


2006 


Jul 13.57 


\^ 1 LUjI nil Uj 


ACIS-S 


(HETG) 




1^1 


2006 


Feb 08.85 


( .n HTt n T'H 

yjl LUjI illl Uj 


ACIS-S 


(HETG) 


RHQ9 


1 5 


2006 


Mar 05.99 


iU/l LUjI Uj 


ACIS-S 


(HETG) 


UOc/O 


±(J.U 


2006 


Apr 01.44 


yjlhUii tUil Uj 


ACIS-S 


(HETG) 


6894 


1 ^ 
xu.u 


2006 


May 14.54 

•J 


\^\vUj\ tUI U 


ACIS-S 


(HETG) 


6896 


1 5 


2006 


Jun 09.76 


f ,n nifi n tti 

yyfbUjI LUjI Uj 


ACIS-S 


(HETG) 


6897 


1 5 


2006 


Aug 12.68 


( n nil n t^n 

1 LUjI IUjI Uj 


ACIS-S 


(HETG) 


RQDI 

Uc/U 1 




2006 


Jun 28.98 


( .n n'Yi n TTI 

K^l (jU/I LUjI Uj 


ACIS-S 


(HETG) 


UOC/O 


14 Q 


2006 


Apr 24.35 


\^l LUjI LUjI U 


ACIS-S 


(HETG) 


68Q5 


14 7 


2006 


Jul 28.46 


( .n n'Yi n Tn 
K^l liU/l LUjI U 


ACIS-S 


(HETG) 


6Q00 


14 6 


2005 


Jul 03.06 


( hi n'Yi n Tn 

KJlbUjllUjI U 


ACIS-S 


(NONE) 




1 9 9 


2005 


Jul 06.33 


1 ,n nn n Tn 
yjivU/l lUi U 


ACIS-S 


(NONE) 




1 9 9 


2005 


Jun 18.48 


( ,n nti n Tn 
yyivU/l vUi U 


ACIS-S 


(NONE) 


5943 


1 9 9 


2005 


Jun 26.90 


( ,n nifi n Tn 
y^ivU/l LUi U 


ACIS-S 


(NONE) 


5946 


1 9 9 


2005 


Jun 01.36 


(In nri nTfi 


ACIS-S 


(NONE) 


5937 


12.2 


2005 


Jun 09.29 


Chandra 


ACIS-S 


(NONE) 


5940 


12.1 


2005 


Jun 15.06 


Chandra 


ACIS-S 


(NONE) 


5942 


12.1 


2005 


Jun 03.95 


Chandra 


ACIS-S 


(NONE) 


5938 


12.0 


2005 


Jun 06.64 


Chandra 


ACIS-S 


(NONE) 


5939 


12.0 


2005 


Jun 11.89 


Chandra 


ACIS-S 


(NONE) 


5941 


12.0 


2005 


Jun 21.22 


Chandra 


ACIS-S 


(NONE) 


5944 


12.0 


2005 


Jun 24.25 


Chandra 


ACIS-S 


(NONE) 


5945 


11.7 


2005 


May 28.83 


Chandra 


ACIS-S 


(NONE) 


5936 


11.6 


2005 


May 26.16 


Chandra 


ACIS-S 


(NONE) 


5935 


11.1 
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Table 1 — Continued 



Date of 


Mission 


Instrument 


Observation 


Exposure 


ol)ser\"atioii 






ID 


(ks) 


2005 Jun 29.56 


Chandra 


ACIS-S (NONE) 


5947 


10.8 


2005 Apr 21.00 


Swift 


XRT 


00035059001 


1.6 


2005 Aug 25.05 


Swift 


XRT 


00035059002 


5.1 


2006 Jun 24.00 


Swift 


XRT 


00035059003 


4.2 


2006 Nov 18.04 


Swift 


XRT 


00035059004 


21.2 


2008 Jan 09.04 


Swift 


XRT 


00036557001 


4.4 


2008 Jan 11.59 


Swift 


XRT 


00036557002 


3.6 


2008 Jan 15.47 


Swift 


XRT 


00036557003 


1.4 


2008 Feb 01.00 


Chandra 


ACIS-S (NONE) 


9122 


10.0 
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Table 2. Unabsorbed X-ray luminosities of SN 1993J at various epochs 



Date of 


Days since 


instrument 




LurnmositA' (10 


erg) 


(il )S('r\'al ioii 






{}.'■-> 2. 1 k('\" 


2 S IvcN" 


(}.'■-> S kr\' 


1993 Apr 03.41 


6.61 


ROSAT-PSPC 


20 63 -1- 1 20 






1993 Apr 05.25-Apr 06.04 


8.45-9.24 


ASCA 


20.50 ± 5.99 


133.02 ± 20.27 


166.85 ± 21.29 


1993 Anr 07 2f)-Ar)r 07 94 


10.45-11.14 


ASCA 


1 9 71 -1- 3 00 


127 48 ± 13 31 


156 91 + 10 72 


1993 Anr OS 34-Ar)r 09 35 


11.54-12.55 


ROSAT-PSPC 


18 76 -1- 92 






1993 Apr 12.23-Apr 13.01 


15.43-16.24 


ROSAT-PSPC 


1 5 63 -1- 76 






1993 Anr 16 83 


20.03 


POSAT-PSPC 


14.96 ± 2.52 






1993 Apr 16.94-Apr 18.73 


20.14-21.93 


ASCA 


13.56 ± 3.15 


69.27 ± 8.32 


87.37 ± 5.36 


1993 Apr 17.68-Apr 19.68 


20.88-22.88 


ROSAT-HRI 


16.54 ± 0.96 






1993 Apr 22.05-Apr 24.0 


25.25-27.20 


ROSAT-PSPC 


14 20 -1- 63 






1993 Apr 25.76 


28.96 


ASCA 




52.91 ± 11.62 


78.22 ± 10.72 


1993 May 01.93-May 02.59 


35.13-35.79 


ASCA 


14.82 ± 8.52 


51.43 ± 8.68 


67.97 ± 14.51 


1993 May 04.30-May 06.30 


37 50-39 50 


ROSAT-PSPC 


1 3 66 -1- 54 






1993 May 12.85-May 13.48 


46.05-46.68 


ROSAT-HRI 


12.47 ± 1.07 






1993 Mav 18 87-Mav 19 86 


52 07-53 06 


ASCA 


1 3 88 -1- 6 15 


29 16±3 97 


50 15 + 14 04 


1993 Oct 24.60-Oct 25.59 


203.30-204.29 


ASCA 




3.84 ± 2.74 




1993 Nov 01.68 


211.38 


ROSAT-PSPC 


7.47 ± 0.23 






1 993 Nov 07 90 


217.60 


ROSAT-PSPC 


8 08 -1- 50 






1994 Apr 01.45 


362.15 


ROSAT-PSPC 


3.59 ± 0.52 






1994 Apr 01.68 


362.38 


ASCA 


8.83 ± 3.15 


2.30 ± 2.06 


11.51 ± 5.36 


1994 Oct 19.17-Oct 21.31 


562.82-564.96 


ROSAT-HRI 


4.35 ± 0.36 






1994 Oct 21.18-Oct 22.33 


564.83-565.98 


ASCA 


4.89 ± 2.37 


< 0.47 


4.89 ± 1.89 


1995 Apr 13.71-May 04.90 


739.36-763.70 


ROSAT-HRI 


5.35 it 0.43 






1995 Oct 19.00 


909.65 


ROSAT-HRI 


6.74 ± 0.48 






1996 Apr 15.88-May 07.83 


1106.53-1144.48 


ROSAT-HRI 


6.59 ± 0.48 






1996 Nov 04.50 


1309.15 


ROSAT-HRI 


5.54 ± 0.80 






1997 Mar 31.10 


1455.75 


ROSAT-HRI 


6.26 it 0.46 






1997 Sep 30.44-Oct 16.00 


1639.09-1654.65 


ROSAT-HRI 


5.83 it 0.43 






1998 Mar 26.10 


1815.75 


ROSAT-HRI 


4.08 it 0.50 






2000 Mar 21.03 


2541.68 


Chandra ACIS-S 


4.44 it 1.40 


1.75 ±0.44 


5.89 ±1.69 


2000 May 07.23 


2588.88 


Chandra ACIS-S 


3.88 ± 0.14 


1.66 ± 0.18 


5.26 ± 0.40 


2001 Apr 22.32 


2938.97 


XMM-Newton EPIC-PN 


4.54 ±0.25 


1.25 ± 0.24 


5.55 ± 0.44 


2005 Feb 24.29-Aug 14.41'» 


4341.94-4513.06 


Chandra ACIS-S (HETG) 


1.77 ±0.18 


0.43 ± 0.08 


2.10 ±0.20 


2005 Apr 21.00-Aug 25.05*' 


4397.65-4523.70 


Swift-XRT 


< 2.24 


< 0.41 


< 2.46 


2005 May 26.16-Jun 15.07"= 


4432.81-4452.72 


Chandra ACIS-S 


2.16 ±0.11 


0.46 ± 0.02 


2.53 ±0.13 


2005 Jun 18.48-Jul 06.33'^ 


4456.13-4473.98 


Chandra ACIS-S 


2.24 ± 0.38 


0.49 ±0.10 


2.63 ± 0.44 


2006 Feb 08.85-Aug 12.68« 


4691.5-4876.33 


Chandra ACIS-S(HETG) 


1.61 ±0.60 


0.24 ±0.07 


1.78 ±0.67 


2006 Jun 24.00-Nov 18. 04^ 


4826.65-4973.69 


Swift-XRT 


1.91 ± 0.32 


0.33 ± 0.05 


2.08 ±0.35 


2008 .Ian 09.04 15.47S 


5385.69-5.392.12 


Sunft-XHT 


2.02 ± 0.55 


0.37 ±0.08 


2.22 ± 0.59 


2008 Fell 01,00 


r)107.()5 


Chaiitlm ACIS-S 


l.()7±().l 


0.:i() ± 0.08 


1.85 = 0.1 



'^Combined observations with IDs 5600, 5601, 6174, 6346 and 6347 
''Combined observations with IDs 00035059001 and 00035059002 
t^Combined observations with IDs 5935, 5936, 5937, 5939, 5940 and 5941 
^Combined observations with IDs 5943, 5944, 5945, 5946, 5947, 5948 and 5949 
''Combined observations with IDs 6892-6901 
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f Combined observations with IDs 00035059003 and 00035059004 

sCombined observations with IDs 00036557001, 00036557002 and 00036557003 



Table 3. Timescales in days for reverse shock to remain radiative for various values of ejecta velocity, mass loss rate, 

composition and density index ^ 



Composition 


Density 
index (n) 


M_5 = 1 


= 0.5 

M-5 =4 


M_5 = 10 


M-5 = l 


= 1-0 

M_5 =4 


M_5 = 10 


M_5 = 1 


Vi = 2.0 
M_5 = 4 


M_5 = 10 


Solar 


7 


60 


250 


570 


4 


16 


40 


0.4 


1.6 


4 




12 


3080 


1.2 X lO'^ 


3 X 10^ 


160 


630 


1580 


8 


32 


80 




20 


1.3 X 10^ 


5.2 X 10^ 


1.3 X 10^ 


3000 


1.2 X 10" 


2.9 X 10" 


125 


500 


1250 


Helium 


7 


190 


760 


1900 


19 


76 


190 


2 


8 


20 




12 


1.2 X lO'^ 


4.9 X 10* 


1.2 X lO-"^ 


473 


1900 


4730 


40 


160 


400 




20 


4.4 X 10^ 


1.8 X 10*^ 


4.4 X 10^ 


1.2 X 10"^ 


4.9 X 10" 


1.2 X 10-'' 


400 


1700 


4000 


Oxygen 


7 


6310 


2.5 X lO'' 


6.3 X 10^ 


470 


1900 


4700 


42 


170 


420 




12 


5.5 X 10^ 


2.2 X 10^ 


5.5 X 10^ 


1.7 X 10" 


6.6 X 10" 


1.7 X 10-^ 


977 


3900 


9770 




20 


3 X 10^ 


1.2 X 10^ 


3 X 10^ 


5.5 X 10^ 


2.2 X lO'^ 


5.5 X 10^ 


1.5 X 10" 


5.9 X 10" 


1.5 X 10^ 



^In this Table 14 is the ejecta velocity in lO" km s ^, and M_5 is mass loss rate in units of 10 ^ Mq yr 
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Table 4. H-a luminosities for SN 1993J 



Days since Luminosity Observatory 
explosion [erg/s] 



Reference 



171 


14.98 X 10^^ 


A m^i fo 


205 


8.67 


X 


1q38 


A m^i fo 


236 


4.96 


X 


1038 


A m^i fo 


255 


4.11 


X 


1038 


A fo 


299 


2.29 


X 


1038 


Asia PT> 


367 


1.28 


X 


1038 


A fo 


553 


2.50 


X 


1038 


Lick 


670 


3.00 


X 


1038 


K^Gck 


687 




X 


iU 


NAOC 


/ UU 


2.75 


X 


1038 


1\T A on 


881 


2.50 


X 


1038 


Lick 


976 


4.80 


X 


1038 


Keck 


986 


2.19 


X 


1038 


NAOC 


998 


1.95 


X 


1038 


NAOC 


1034 


1.99 


X 


1038 


NAOC 


1280 


1.74 


X 


1038 


NAOC 


1318 


1.29 


X 


1038 


NAOC 


1395 


1.38 


X 


1038 


NAOC 


1729 


1.10 


X 


1038 


NAOC 


1766 


0.86 


X 


1038 


Keck 


2066 


0.65 


X 


1038 


NAOC 


2454 


0.37 


X 


1038 


Keck 


3149 


0.50 


X 


1038 


NAOC 


3401 


0.43 


X 


1038 


NAOC 


3503 


0.43 


X 


1038 


NAOC 


3610 


0.28 


X 


1038 


Keck 



Patat et al (1995) 
Patat et al (1995) 
Patat et al (1995) 
Patat et al (1995) 
Patat et al (1995) 
Patat et al (1995) 
Matheson et al. (2000b) 
Matheson et al. (2000b) 
Zhang et al. (2004) 
Zhang et al. (2004) 
Matheson et al. f 2000b) 
Matheson et al. (2000a) 
Zhang et al. (2004) 
Zhang et al. (2004 ) 
Zhang et al. (2004) 
Zhang et al. (2004) 
Zhang et al. (2004) 
Zhang et al. (2004) 
Zhang et al. (2004) 
Matheson et al. (2000b) 
Zhang et al. (2004) 
Matheson et al. (2000b) 
Zhang et al. (2004) 
Zhang et al. (2004) 
Zhang et al. (2004) 
Filippenko k Matheson (2005) 
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Fig. 1. — Here we plot the two component Mekai fit to SN 1993J Chandra data of 2005 taken 
with ACIS-S. Best fit temperatures from the Mekal model are 0.73 keV and 2.21 keV with 
absorption column density of 6 x 10^ 



cm-2. 



-so- 




lo 100 1000 10 100 1000 



Days since explosion Days since explosion 




Days since explosion Days since explosion 

Fig. 2. — The X-ray light curves of SN 1993 J in the 0.3-8 keV (upper left panel), 2- 
8 keV (upper right panel), and 0.3-2.4 keV (lower left panel) bands observed with multiple 
telescopes. For around first 1000 days, 0.3-2.4 keV light curves declines slowly (t^°'^^) 
whereas the 2-8 keV light curves decline as t~^. The overall 0.3-8 keV light curve declines 
as indicating radiative nature of the reverse shock. After around day 1000, the shocks 
seems to become adiabatic. The lower right panel shows the comparison between the 0.3- 
2.4 keV and 2-8 keV components and demonstrates that the soft component take overs 
around day 200 and dominates at late epochs. 
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SN 1993J Hard X-ray Luminosity 



Data Points 
Hydro Calculation 




10 100 1000 10000 

Time (Days) 

Fig. 3. — Model plot of the Hard X-ray band (1-6 A) luminosity evolution calculated 
from liydrodynamic modelling with VH-1 and CHIANTI codes starting from a Nomoto 
and Shigeyama 4H47 model. We plot this model with the hard band observed luminosities 
in 2-8 keV band. The model does not trace the early luminosity evolution for first 10-15 
days (probably due to complicated X-ray emission process at early phase), however seems 
to reproduce the observed light curve well after that. 
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Q Q I I I I I ■ I I I I I I I ■ I I I I Li I I I I ■ I I I I I I i_ 

10 100 1000 

Days since explosion 

Fig. 4. — Here we plot the ratio of 2-8 keV luminosities versus 0.3-2.4 keV luminosities 
at various epochs. The hardness ratio decreases with time and roughly follows a powerlaw 
dependence with time index of —0.55. 
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H-alpha plotted with 2.0-8.0 keV light curve 
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Fig. 5. — Plot of Ha luminosity light curve (taken from Table HI). We also plot the 2-8 keV 
X-ray light curve. The Ha roughly traces the 2-8 keV band light curve after day ~ 300 
indicating that Ha arises due to the CS interaction. The Ha luminosities are high fractions 
(30-50%) of X-ray luminosities, which probably indicate towards the presence of clumps in 
the ejecta. 



